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Abstract. A new analytical approach to the motion 
and radiation of (comparable mass) binary systems has 
been introduced in 1999 under the name of Effective 
One Body (EOB) formalism. We review the basic 
elements of this formalism, and discuss some of its recent 
developments. Several recent comparisons between EOB 
predictions and Numerical Relativity (NR) simulations have 
shown the aptitude of the EOB formalism to provide accurate 
descriptions of the dynamics and radiation of various binary 
systems (comprising black holes or neutron stars) in regimes 
that are inaccessible to other analytical approaches (such as 
the last orbits and the merger of comparable mass black holes). 
In synergy with NR simulations, post-Newtonian (PN) theory 
and Gravitational Self-Force (GSF) computations, the EOB 
formalism is likely to provide an efficient way of computing 
the very many accurate template waveforms that are needed 
for Gravitational Wave (GW) data analysis purposes. 



1. Introduction 

The general relativistic iV-body problem has been investigated from the 
early days of Einstein's gravitation theory (and even earlier, because it 
was already tackled by Johannes Droste within the framework of the 
1913 Einstein- Grossmann "Entwurf" theory). Here, we shall focus on 
the general relativistic two-body problem. This problem has been the 
subject of many investigations within the post- Newtonian (PN) formalism, 



since the pioneering works of Einstein (1915; when mi <C TO2), Lorentz 
and Droste (1917), Levi-Civita (1937) and Einstein, Infeld and Hoffmann 
(1938). [see, e.g., [IS] for a review and references to the early literature.] 
For many years, the first post- Newtonian (1PN) approximation (i.e. 
the inclusion of the leading-order relativistic corrections, proportional to 
(v/c) 2 or GM/c 2 r, to the Newtonian equations of motion) appeared as 
being accurate enough for applying Einstein's theory to known binary 
systems. The situation changed in the mid 1970's with the discovery of 
the Hulse- Taylor binary pulsar PSR 1913 + 16. The need to compare 
the accurate observations of this system (by Taylor and collaborators) 
to the predictions of Einstein's theory motivated the development of 
improved relativistic theories of binary systems, applicable to strongly self- 
gravitating bodies, and including terms up to the 2.5PN approximation 
(i.e. 0[(v/c) 5 ) beyond Newton). [See [H] and references therein.] 
The situation has again changed recently with the development of 
interferometric gravitational wave (GW) detectors, and the prospect of 
detecting the GW's emitted during the last orbits and the coalescence of 
binary systems made of black holes or neutron stars. The latter prospect 
motivated the development (or improvement) of several different methods 
of computing the motion and radiation of binary systems. 

First, this motivated pushing PN calculations of the dynamics of 
binary systems to the 3PN level [U] GUI Ml [Ml W\, with inclusion of 
3.5PN radiation- reaction terms [1051 [531 155] . Second, this motivated 
the development of new, accurate GW generation formalisms, notably 
the Blanchet-Damour-Iyer (matched) "multipolar post-Minkowskian" 
formalism [211 G2] [55] [56] [23] [17] and the "direct integration of the relaxed 
Einstein's equations" formalism of Will and collaborators |119l 11181 1104J . 
which extended previous work by Epstein and Wagoner [78] and Thorne 
[116J - These GW generation formalisms allowed one to compute emitted 
gravitational waveforms with an unprecedented PN accuracjo After the 
1PN correction to the waveform [78] Ell [22] , there is a 1.5PN "tail" (i.e. 
hereditary) correction [231 USD] EH] , then a "direct" 2PN term [28112711115] . 
followed by higher-order corrections [IH][3J[3J[2J[251[in][Ill[Sl[3T]. 
[See [2D] for a detailed account and more references.] Parallely to these 
improved PN computations of the GW emission of comparable-mass 
systems (with mi ~ m 2 ), other authors developed the analytical theory 
of the GW emission of extreme mass-ratio systems (with mi <C 1712): see 
Refs. [Ml rnUJ [333] HH] and the review of Sasaki and Tagoshi [TIT] . 

1 For gravitational waveforms, one conventionally defines the PN accuracy as the 
fractional PN accuracy with respect to the leading-order, 0(c ), quadrupolar 
emission. E.g., a lPN-accurate waveform retains next-to- leading order terms, 
i.e. terms smaller than the quadrupolar waveform by a factor O(o~ 2 ). 



Some of the PN calculations of the dynamics, and/or GW emission, 
of comparable-mass systems have been recently (re)done (e.g. the 3PN 
dynamics 80 ) by using a somewhat different formalism, dubbed "effective 
field theory" [85]. Let us, however, note that most of the technical aspects 
of the effective- field-theory approach had already been introduced and used 
before. For instance: (i) Ref. [50 discussed the (Fokker) two-body effective 
action due to the exchange of a linear field (of spin s = 0, 1 and 2) ; 
(ii) Ref. |51j explicitly discussed the representation (and computation) 
of the (nonlinear) effective two-body action in terms of Feynman-like 
diagrams (made of concatenated propagators and vertices); (iii) The 
appendix A of Ref. 52 discussed finite-size effects in terms of nonminimal 
worldline couplings in the effective action; (iv) The (quantum field theory) 
technique of dimensional regularization (together with a diagrammatic 
analysis of ultraviolet divergences) was crucially used to derive the 3PN 
dynamics in Refs. (62J [24], and 3PN radiation in Ref. [26]; and (v) The 
exponential parametrization of the metric (which suppresses the leading- 
order gravitational cubic vertex) had been introduced in Ref. [35] and 
then standardly used in many PN works. It is, however, possible that the 
more systematic (and automatized) implementation of such diagrammatic 
methods, together with the tapping of standard techniques for computing 
Feynman graphs (as exemplified in |80| ) may allow one to be more efficient 
in computing higher-order processes, or, at least, to open new ways of 
understanding them (see, in this respect, Ref. [HI])- 

Separately from these purely analytical approaches to the motion and 
radiation of binary systems, which have been developed since the early days 
of Einstein's theory, Numerical Relativity (NR) simulations of Einstein's 
equations have relatively recently (2005) succeeded (after more than thirty 
years of developmental progress) to stably evolve binary systems made of 
comparable mass black holes |107| l44ll4"ll35 j . This has led to an explosion of 
works exploring many different aspects of strong-field dynamics in General 
Relativity, such as spin effects, recoil, relaxation of the deformed horizon 
formed during the coalescence of two black holes to a stationary Kerr black 
hole, high- velocity encounters, etc.; see 108 for a review. In addition, 
recently developed codes now allow one to accurately study the orbital 
dynamics, and the coalescence of binary neutron stars. Much physics 
remains to be explored in these systems, especially during and after the 
merger of the neutron stars (which involves a much more complex physics 
than the pure-gravity merger of two black holes). 

Recently, a new source of information on the general relativistic two- 
body problem has opened: gravitational self-force (GSF) theory. This 
approach goes one step beyond the test-particle approximation (already 
used by Einstein in 1915) by taking into account self-field effects that 



modify the leading-order geodetic motion of a small mass m\ moving in the 
background geometry generated by a large mass ni2. After some ground 
work (notably by DeWitt and Brehme) in the 1960's, GSF theory has 
recently undergone rapid developments (mixing theoretical and numerical 
methods) and can now yield numerical results that yield access to new 
information on strong-field dynamics in the extreme mass-ratio limit 
Tfii rri2. See Ref. [3] for a review. 

Each of the approaches to the two-body problem mentioned so far, 
PN theory, NR simulations and GSF theory, have their advantages and 
their drawbacks. It has become recently clear that the best way to 
meet the challenge of accurately computing the gravitational waveforms 
(depending on several continuous parameters) that are needed for a 
successful detection and data analysis of GW signals in the upcoming 
LIGO/Virgo/GEO/. . . network of GW detectors is to combine knowledge 
from all the available approximation methods: PN, NR and GSF. Several 
ways of doing so are a priori possible. For instance, one could try to 
directly combine PN-computed waveforms (approximately valid for large 
enough separations, say r > 10G(mi + m2)/c 2 ) with NR waveforms 
(computed with initial separations r$ > 10G(mi + W2)/c 2 and evolved 
up to merger and ringdown). However, this method still requires too 
much computational time, and is likely to lead to waveforms of rather 
poor accuracy, see, e.g., [551 1571 l§7] . 

On the other hand, five years before NR succeeded in simulating the 
late inspiral and the coalescence of binary black holes, a new approach 
to the two-body problem was proposed: the Effective One Body (EOB) 
formalism [JUl SU ED II7J- The basic aim of the EOB formalism is to 
provide an analytical description of both the motion and the radiation of 
coalescing binary systems over the entire merger process, from the early 
inspiral, right through the plunge, merger and final ringdown. As early as 
2000 [41] this method made several quantitative and qualitative predictions 
concerning the dynamics of the coalescence, and the corresponding GW 
radiation, notably: (i) a blurred transition from inspiral to a 'plunge' that 
is just a smooth continuation of the inspiral, (ii) a sharp transition, around 
the merger of the black holes, between a continued inspiral and a ring-down 
signal, and (iii) estimates of the radiated energy and of the spin of the 
final black hole. In addition, the effects of the individual spins of the black 
holes were investigated within the EOB [ITJ ISFJ and were shown to lead to 
a larger energy release for spins parallel to the orbital angular momentum, 
and to a dimensionless rotation parameter J/E 2 always smaller than 
unity at the end of the inspiral (so that a Kerr black hole can form 
right after the inspiral phase). All those predictions have been broadly 
confirmed by the results of the recent numerical simulations performed 



by several independent groups (for a review of numerical relativity results 
and references see |108j ). Note that, in spite of the high computer power 
used in NR simulations, the calculation, checking and processing of one 
sufficiently long waveform (corresponding to specific values of the many 
continuous parameters describing the two arbitrary masses, the initial spin 
vectors, and other initial data) takes on the order of one month. This is a 
very strong argument for developing analytical models of waveforms. 

2. EOB description of the conservative dynamics of two body 



Before reviewing some of the technical aspects of the EOB method, let us 
indicate the historical roots of this method. First, we note that the EOB 
approach comprises three, rather separate, ingredients: 

(i) a description of the conservative (Hamiltonian) part of the dynamics 
of two bodies; 

(ii) an expression for the radiation-reaction part of the dynamics; 

(iii) a description of the GW waveform emitted by a coalescing binary 



For each one of these ingredients, the essential inputs that are used in 
EOB works are high-order post-Newtonian (PN) expanded results which 
have been obtained by many years of work, by many researchers (see 
references above). However, one of the key ideas in the EOB philosophy 
is to avoid using PN results in their original "Taylor-expanded" form (i.e. 
c + ci v/c+c 2 v 2 /c 2 + c 3 v 3 /c 3 + - ■■+c n v n /c n ), but to use them instead in 
some resummed form (i.e. some non-polynomial function of v/c, defined 
so as to incorporate some of the expected non-perturbative features of 
the exact result). The basic ideas and techniques for resumming each 
ingredient of the EOB are different and have different historical roots. 

Concerning the first ingredient, i.e. the EOB Hamiltonian, it was 
inspired by an approach to electromagnetically interacting quantum two- 
body systems introduced by Brezin, Itzykson and Zinn- Justin [36 . 




The resummation of the second ingredient, i.e. the EOB radiation- 
reaction force J 7 , was initially inspired by the Pade resummation of 
the flux function introduced by Damour, Iyer and Sathyaprakash [58] . 
More recently, a new and more sophisticated resummation technique 
for the radiation reaction force J- has been introduced by Damour and 
Nagar [57l [67]. 

As for the third ingredient, i.e. the EOB description of the waveform 
emitted by a coalescing black hole binary, it was mainly inspired by the 
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work of Davis, Ruffini and Tiomno [77] which discovered the transition 
between the plunge signal and a ringing tail when a particle falls into a 
black hole. [Additional motivation for the EOB treatment of the transition 
from plunge to ring-down came from work on the, so-called, "close limit 
approximation" [109] .] In addition, a very efficient resummation of the 
waveform has been introduced by Damour, Iyer and Nagar 64, 66, 57 . It 
will be discussed in detail below. 

Within the usual PN formalism, the conservative dynamics of a two- 
body system is currently fully known up to the 3PN level |90l |30l |62j 
[24j l89l [80] (see below for the partial knowledge beyond the 3PN level) . 
Going to the center of mass of the system (p 1 + p 2 =0), the 3PN-accurate 
Hamiltonian (in Arnowitt-Deser-Misner-type coordinates) describing the 
relative motion, q = q\ — <72, p = pi = —P2, has the structure 

H$f ve (q, P ) = H {q,p) + Mq,p) + ^ h^p) + ^ Mq.p) , (i) 

where 
with 

M = mi + ni2 and fx = mi m^jM , (3) 

corresponds to the Newtonian approximation to the relative motion, while 
H2 describes IPN corrections, H4 2PN ones and H§ 3PN ones. In 
terms of the rescaled variables q' = q/GM, p' = p/ fi, the explicit form 
(after dropping the primes for readability) of the 3PN-accurate rescaled 
Hamiltonian H = H/fi reads [501 ESI E2] 

H N {q,p) = ^--- (4) 
z q 



ffiPN(g.p) = 1)(P 2 ) 2 - ^[(3 + v)v 2 + v{n ■ P?]- q + ^ , (5) 
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8 q 

+ i[(5 + 8v)p 2 + 3 V (n ■ p) 2 }± - 1(1 + 3i/)^ , (6) 
H 3PN (q,p) = T ^(-5 + 35 i .-70^ + 35^ 3 )(p 2 ) 4 
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In these formulas v denotes the symmetric mass ratio: 
_ [i _ m\ TO2 



M (mi+m 2 ) 



2 ' 



(7) 



(8) 



The dimensionless parameter v varies between (extreme mass ratio case) 
and \ (equal mass case) and plays the role of a deformation parameter 
away from the test-mass limit. 

It is well known that, at the Newtonian approximation, Ho(q,p) can 
be thought of as describing a 'test particle' of mass fi orbiting around 
an 'external mass' GM. The EOB approach is a general relativistic 
generalization of this fact. It consists in looking for an 'effective external 
spacetime geometry' g c ^(x x \GM,v) such that the geodesic dynamics of 
a 'test particle' of mass /i within g c ^{x x ,GM,v) is equivalent (when 
expanded in powers of 1/c 2 ) to the original, relative PN-expanded 
dynamics |T]). 

Let us explain the idea, proposed in jJD], for establishing a 'dictionary' 
between the real relative- motion dynamics, ([lj, and the dynamics of 
an 'effective' particle of mass /i moving in jf^ (x x , GM, v) . The idea 
consists in 'thinking quantum mechanicallyQ Instead of thinking in 
terms of a classical Hamiltonian, H(q,p) (such as H^p^ lve , Eq. (p}), 
and of its classical bound orbits, we can think in terms of the quantized 



2 This is related to an idea emphasized many times by John Archibald Wheeler: 
quantum mechanics can often help us in going to the essence of classical 
mechanics. 



energy levels E(n,£) of the quantum bound states of the Hamiltonian 
operator H(q,p). These energy levels will depend on two (integer valued) 
quantum numbers n and I. Here (for a spherically symmetric interaction, 
as appropriate to fj lelatlve ^ g parametrizes the total orbital angular 
momentum (L 2 — £[£ + 1) h 2 ), while n represents the 'principal quantum 
number' n = £+n r +l, where n r (the 'radial quantum number') denotes the 
number of nodes in the radial wave function. The third 'magnetic quantum 
number' m (with —£ < m < £) does not enter the energy levels because of 
the spherical symmetry of the two-body interaction (in the center of mass 
frame). For instance, the non-relativistic Newton interaction Eq. |2]) gives 
rise to the well-known result 

which depends only on n (this is the famous Coulomb degeneracy). When 
considering the PN corrections to Ho, as in Eq. ([lj, one gets a more 
complicated expression of the form 



r-i< Mauve / /n 

£ 3 pn (n,£) = -~^ 



— (211 + C20 \ 
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(10) 



where we have set a = GM/i/H = G mi 1112 /h, and where we consider, 
for simplicity, the (quasi-classical) limit where n and £ are large numbers. 
The 2PN-accurate version of Eq. (fTU]) had been derived by Damour and 
Schafer [76] as early as 1988 while its 3PN-accurate version was derived 
by Damour, Jaranowski and Schafer in 1999 [60]. The dimensionless 
coefficients c pq are functions of the symmetric mass ratio v = /i/M, for 
instance C40 = ^ (145 — \hv + v 2 ). In classical mechanics (i.e. for large 
n and £), it is called the 'Delaunay Hamiltonian', i.e. the Hamiltonian 
expressed in terms of the action variable^ J = £h = § p v dip, and 
N = nh = I r + J, with I r = y~ § p r dr. 

The energy levels (|10p encode, in a gauge-invariant way, the 3PN- 
accurate relative dynamics of a 'real' binary. Let us now consider an 
auxiliary problem: the 'effective' dynamics of one body, of mass \x, 
following (modulo the Q term discussed below) a geodesic in some in- 
dependent 'effective external' (spherically symmetric) metric0 

gf v dx^ dx v = -A(R;v)c 2 dT 2 + B{R;v)dR 2 + R 2 {d6 2 + sin 2 9 dip 2 ) . (11) 

3 We consider, for simplicity, 'equatorial' motions with m — £ , i.e., classically, 
= f. 

4 It is convenient to write the 'effective metric' in Schwarzschild-like coordinates. 
Note that the effective radial coordinate R differs from the two-body ADM- 



Here, the a priori unknown metric functions A(R; v) and B(R; v) will be 
constructed in the form of expansions in GM/c 2 R: 



„,„ v , ~ GM „ / GM \ 2 „ /GM\ 3 _ / GAf \ 4 

r GAf ~ / GAf \ 2 , fGM\ 3 , . 

B(R;v) = i + 6l _ +6a ^_j ^(—J +..-, (12) 

where the dimensionless coefheients a n ,b n depend on From the 
Newtonian limit, it is clear that we should set a\ = —2. In addition, 
as v can be viewed as a deformation parameter away from the test-mass 
limit, we require that the effective metric tend to the Schwarzschild 
metric (of mass M) as v — > 0, i.e. that 

A(R; v = 0) = 1 - 2GM/c 2 R = B^ 1 (R\ v = 0) . 



Let us now require that the dynamics of the "one body" [i within the 
effective metric <?°^ be described by an "effective" mass-shell condition of 
the form 



n^ v n cS r) cff 
dcfiPfi Pv 



1 2 + Q(pf) = o, 



where Q{p) is (at least) quartic in p. Then by solving (by separation of 
variables) the corresponding 'effective' Hamilton- Jacobi equation 
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(13) 



one can straightforwardly compute (in the quasi-classical, large quantum 
numbers limit) the effective Delaunay Hamiltonian £ e g (N e g , J e s), with 
Ncs = n cS h, Jeff = 4ff h (where iV off = J cff + if, with lf = ^§ pf dR, 
dS e g(R)/dR). This yields a result of the form 
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(14) 



coordinate relative distance R = |qf|. The transformation between the two 
coordinate systems has been determined in Refs. [401 161] . 



where the dimensionless coefficients are now functions of the unknown 
coefficients a n ,b n entering the looked for 'external' metric coefficients (|12l) . 

At this stage, one needs to define a 'dictionary' between the real 
(relative) two-body dynamics, summarized in Eq. fjlO[) . and the effective 
one-body one, summarized in Eq. (| 14[) . As, on both sides, quantum 
mechanics tells us that the action variables are quantized in integers 
(■Wreai = N c ff = n D f[h, etc.) it is most natural to identify n = n e s 
and £ = £ c g. One then still needs a rule for relating the two different 
energies E^ tlvc and £ e g. Ref. [40] proposed to look for a general map 
between the real energy levels and the effective ones (which, as seen when 
comparing ([TU)) and (|14l) , cannot be directly identified because they do not 
include the same rest- mass contribution^, namely 



c „ / ^relative \ ^relative / ^relative / ^relative 
tcS 1 f fc rcal \ A'cal I , „ fc rcal , „ I fc rcal 
2 ~ 1 = J \ 2 = 2 1 + "1 5 h a2 2 



(^relative \ ^ \ 
-ft-) + -)- <15) 

The 'correspondence' between the real and effective energy levels is 
illustrated in Fig. [TJ 

Finally, identifying £ e $ (n, £)/fic 2 to 1 + f{E^ tivc (n,£)/fic 2 ) yields 
a system of equations for determining the unknown EOB coefficients 
a n ,b n ,a n , as well as the three coefficients zi, Za, z% parametrizing a general 
3PN-level quartic mass-shell deformation: 

<?3Pn(p) = (^jf) [zip 4 + z 2 p 2 (n-p) 2 + z 3 (n-p) 4 ] . 

[The need for introducing a quartic mass-shell deformation Q only arises 
at the 3PN level.] 

The above system of equations for a n ,b n , a n (and z% at 3PN) was 
studied at the 2PN level in Ref. [40], and at the 3PN level in Ref. [ST] . 
At the 2PN level it was found that, if one further imposes the natural 
condition b\ = +2 (so that the linearized effective metric coincides with 
the linearized Schwarzschild metric with mass M = m\ +W2), there exists 
a unique solution for the remaining five unknown coefficients 02, 03, 62, «i 
and «2- This solution is very simple: 

O2=0, a 3 = 2i/, b 2 = 4-6z/, ax = ^ , a 2 = 0. (16) 

5 Indeed E%£f = Mc 2 + E T ^ tivc = Mc 2 + Newtonian terms + lPN/c 2 + ■ ■ ■, 
while £" cffcctivc = fic 2 + N + lPN/c 2 + • ■ •. 



e = f(E) 




Figure 1. Sketch of the correspondence between the quantized energy 
levels of the real and effective conservative dynamics, n denotes the 
'principal quantum number' (n = n r + £ + 1, with n r = 0, 1, . . . denoting 
the number of nodes in the radial function), while I denotes the (relative) 
orbital angular momentum (L 2 = £(£ + 1) h 2 ). Though the EOB method is 
purely classical, it is conceptually useful to think in terms of the underlying 
(Bohr-Sommerfeld) quantization conditions of the action variables Ir and 
J to motivate the identification between n and £ in the two dynamics. 



At the 3PN level, it was found that the system of equations is consistent, 
and underdetermined in that the general solution can be parametrized by 
the arbitrary values of z\ and z%. It was then argued that it is natural to 
impose the simplifying requirements z\ — = zi , so that Q is proportional 
to the fourth power of the (effective) radial momentum p r . With these 
conditions, the solution is unique at the 3PN level, and is still remarkably 
simple, namely 

04 = a4 v , ds = 2(3^ — 26)^ , = , z 3 = 2(4 — Zv)v . 

Here, denotes the number 

94 41 

a 4 = 7T 2 ~ 18.6879027 , (17) 

3 32 y J 

while g?3 denotes the coefficient of (GM/c 2 R) 3 in the PN expansion of the 
combined metric coefficient 

D(R) = A(R) B(R) . 



Replacing B(R) by D(R) is convenient because (as was mentioned above), 
in the test-mass limit v — > 0, the effective metric must reduce to the 



Schwarzschild metric, namely 



A(R;is = 0) = B- 1 (R;v = 0) = 1-2^0 , 

so that 

D(R;v = 0) =1. 

The final result is that the three EOB potentials A, D, Q describing the 
3PN two-body dynamics are given by the following very simple results. In 
terms of the EOB "gravitational potential" 

GM 



~ c 2 R ' 

A 3PN (R) = I - 2u + 2vu 3 + a i vu 4 , (18) 
D 3PN (R) = {A(R)B(R)) 3PN = 1 - Qvu 2 + 2(3// - 26)//u 3 , (19) 



Q 3 Mq lP ) = ^2{A-Z V )v u 2 4- (20) 

In addition, the map between the (real) center-of-mass energy of the 
binary system E™^ ive = # relative = £^ tive - M c 2 and the effective one 
£ e ff is found to have the very simple (but non trivial) form 

„ ^relative / ,, ^relative \ „ ™2 „4 ™2 „4 

% = l + ^-(l + ?%-)= oT nT^ (21) 



where s = {£ * c °^) 2 = (Afc 2 + ££*f ive ) 2 is Mandelstam's invariant s = 
-(Pi +P2) 2 . 

It is truly remarkable that the EOB formalism succeeds in condensing 
the complicated, original 3PN Hamiltonian, Eqs. (d))-©, into the 
very simple potentials A, D and Q displayed above, together with the 
simple energy map Eq. (|21jl. For instance, at the 1PN level, the 
already somewhat involved Lorentz-Droste-Einstein-Infeld-Hoffmann 1PN 
dynamics (Eqs. (jlj and (|5]l) is simply described, within the EOB 
formalism, as a test particle of mass fi moving in an external Schwarzschild 
background of mass M = mi + m2, together with the (crucial but quite 
simple) energy transformation (|21[) . [Indeed, the ^-dependent corrections 
to A and D start only at the 2PN level] At the 2PN level, the seven 
rather complicated //-dependent coefficients of H 2P N{q,p), Eq. ([5]), get 
condensed into the two very simple additional contributions + 2z/u 3 in 
A(u), and — Qvu 2 in D(u). At the 3PN level, the eleven quite complicated 
independent coefficients of H 3P ^, Eq. ([7]), get condensed into only three 



simple contributions: + a^vu 4 in A(u), + 2(3z/ — 26)vu 3 in D(u), and Q3PN 
given by Eq. (|20p . This simplicity of the EOB results is not only due 
to the reformulation of the PN-expandcd Hamiltonian into an effective 
dynamics. Indeed, the A-potential happens to be much simpler that 
it could a priori have been: (i) as already noted it is not modified at 
the IPN level, while one would a priori expect to have found a 1PN 
potential -Aipn(^) = 1 — 2u + va,2U 2 with some non zero 02; and (ii) 
there are striking cancellations taking place in the calculation of the 2PN 
and 3PN coefficients 02 (^) and 03(^)5 which were a priori of the form 
a,2 {v) = 0,2V + a' 2 v 2 , and 0,3(1/) = a^v + alji/ 2 + a 3 V 3 , but for which the 
^-nonlinear contributions a' 2 i/ 2 ,a' 3 v 2 and a 3 V 3 precisely cancelled out. 

The fact that the 3PN coefficient 04 in the crucial 'effective radial 
potential' ^3pn(-R), Eq. ([15]). is rather large and positive indicates that 
the ^-dependent nonlinear gravitational effects lead, for comparable masses 
(v ~ j), to a last stable (circular) orbit (LSO) which has a higher frequency 
and a larger binding energy than what a naive scaling from the test-particle 
limit (v — > 0) would suggest. Actually, the PN-expanded form (TT5|) of 
^3Pn(-R) does not seem to be a good representation of the (unknown) 
exact function Aeob(^) when the (Schwarzschild-like) relative coordinate 
R becomes smaller than about 6GM/c 2 (which is the radius of the LSO 
in the test-mass limit). By continuity with the test-mass case, one a 
priori expects that A^p^(R) always exhibits a simple zero defining an 
EOB "effective horizon" that is smoothly connected to the Schwarzschild 
event horizon at R — 2GM/c 2 when v — > 0. However, the large value 
of the 04 coefficient does actually prevent ^3pn to have this property 
when v is too large, and in particular when v = 1/4. It was therefore 
suggested [51] to further resunJ! ^4 3 pn(-R) by replacing it by a suitable 
Pade (P) approximant. For instance, the replacement of ^3pn(-R) bj0 

A\ (R) P3 1 [A 3PN (R)]= VT 1 "^ 3 ( 22 ) 

ensures that the v = -7 case is smoothly connected with the v = limit. 

The same kind of ^-continuity argument, discussed so far for the A 
function, needs to be applied also to the Z?3pn(-R) function defined in 

6 The PN-expanded EOB building blocks A 3PN (-R), B :iPN (R), . . . already 
represent a resummation of the PN dynamics in the sense that they have 
"condensed" the many terms of the original PN-expanded Hamiltonian within 
a very concise format. But one should not refrain to further resum the EOB 
building blocks themselves, if this is physically motivated. 

7 We recall that the coefficients m and (di, c/2, ofe) of the (1, 3) Pade approximant 
P3 1 [^3pn(w)] are determined by the condition that the first four terms of the 
Taylor expansion of A\ in powers of u — GM/(c 2 R) coincide with ^pn. 



Eq. (|19p . A straightforward way to ensure that the D function stays 
positive when R decreases (since it is D = 1 when v — > 0) is to replace 
D 3PN {R) by D%(R) = P 3 ° [D 3PN (i?)], where P§ indicates the (0,3) Pade 
approximant and explicitly reads 

D 3(«) = m 2 ^ ( 23 ) 

1 + wu 1 — 2[6v — 2Q)vu i 

3. EOB description of radiation reaction and of the emitted 
waveform during inspiral 

In the previous Section we have described how the EOB method encodes 
the conservative part of the relative orbital dynamics into the dynamics 
of an 'effective' particle. Let us now briefly discuss how to complete 
the EOB dynamics by defining some resummed expressions describing 
radiation reaction effects, and the corresponding waveform emitted at 
infinity. One is interested in circularized binaries, which have lost their 
initial eccentricity under the influence of radiation reaction. For such 
systems, it is enough (in first approximation |41j : see, however, the 
recent results of Bini and Damour 15 ) to include a radiation reaction 
force in the P v equation of motion only. More precisely, we are using 
phase space variables R, Pr, <p, P v associated to polar coordinates (in 
the equatorial plane 9 — -|). Actually it is convenient to replace the 
radial momentum Pr by the momentum conjugate to the 'tortoise' radial 
coordinate R* = J dR(B / A) 1 / 2 7 i.e. P Rt = (A/B)^ 2 P R . The real EOB 
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Hamiltonian is obtained by first solving Eq. (f2~Tj) to get if' 
terms of £ e ff, and then by solving the effective Hamilton- Jacobi equation 
to get £ G ff in terms of the effective phase space coordinates q c g and p c s- 
The result is given by two nested square roots (we henceforth set c = 1): 

rrrcal i i 

H E OB(r,p r „<p) = = -VI + 2v{H eS - 1) , (24) 



where 




(25) 



with z 3 = 2v (4 — 3^). Here, we are using suitably rescaled dimensionless 
(effective) variables: r = R/GM, p r „ = Pr,//U, p v — P^/fiGM, as well 
as a rescaled time t = T/GM. This leads to equations of motion for 
(r, ip,p r ,,Pip) of the form 

% - d -p^ = n, (26) 

dt dp v 



dr fA\ 1/2 8H EO b 



dt \B J dp Tt ' 

dp v 



(27) 



dt \BJ dr 

which explicitly read 

dip _ Ap v 



T v , (28) 

(29) 
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dt vr 2 HH cK 
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n , (30) 



dt \bj „im eff v^ + * 3 ^ J ' (31) 



dp v 
~dT 



K , (32) 



dpr, ( A\ 1/2 1 



dt \BJ 2vHH eS 

{t^-^y^-^yy ,33) 

where A 1 = dA/dr. As explained above the EOB metric function A(r) 
is defined by Pade resumming the Taylor-expanded result (|T2"j) obtained 
from the matching between the real and effective energy levels (as we were 
mentioning, one uses a similar Pade resumming for D{r) = A(r)B(r)). 
One similarly needs to resum J- v , i.e., the tp component of the radiation 
reaction which has been introduced on the r.h.s. of Eq. 



Several methods have been tried during the development of the EOB 
formalism to resum the radiation reaction J- v (starting from the high- 
order PN-expanded results that have been obtained in the literature; see 
references in the Introduction above). Here, we shall briefly explain the 
new, parameter-free resummation technique for the multipolar waveform 
(and thus for the energy flux) introduced in Ref. [53J [S3] and perfected 
in [S7]. To be precise, the new results discussed in Ref. [57] are twofold: 
on the one hand, that work generalized the I = m = 2 resummed 
factorized waveform of |64[ 166] to higher multipoles by using the most 
accurate currently known PN-expanded results [92] [14] [31] as well as the 
higher PN terms which are known in the test- mass limit (1131 1114] ; on 
the other hand, it introduced a further resummation procedure which 
consists in considering a new theoretical quantity, denoted as pi m (%), 
which enters the (£, m) waveform (together with other building blocks, see 
below) only through its £-th power: hg m oc (pi m (x)) e . Here, and below, 



x denotes the invariant PN-ordering parameter given during inspiral by 
x= (GMn/c 3 ) 2 / 3 . 

The main novelty introduced by Refs. [64j [66] [57] is to write the (£, m) 
multipolar waveform emitted by a circular nonspinning compact binary as 
the product of several factors, namely 

h& = ^n&^x^V^ (f , *) Sip im e^4 m - (34) 

Here e denotes the parity of £+m (e = 7r(£+rn)), i.e. e = for "even-parity" 
(mass-generated) multipoles (£ + m even), and e = 1 for "odd-parity" 
(current-generated) ones (£ + m odd); n^L and Q+ e (^) are numerical 

coefficients; S^g is a /^-normalized effective source (whose definition comes 
from the EOB formalism) ; Tg m is a resummed version |64l 166] of an infinite 
number of "leading logarithms" entering the tail effects [23] [18] ; S( m is a 
supplementary phase (which corrects the phase effects not included in the 
complex tail factor Ti m ), and, finally, (pt m ) denotes the £-ih power of the 
quantity pi m which is the new building block introduced in [ST] . Note that 
in previous papers [531 1HS] the quantity (pe m ) £ was denoted as ft m and we 
will often use this notation below. Before introducing explicitly the various 
elements entering the waveform (|34l) it is convenient to decompose h(, m as 

where is the Newtonian contribution (i.e. the product of the first 

five factors in Eq. (|34]l ) and 

= S^T lm e l5 ^f tm (36) 

represents a resummed version of all the PN corrections. The PN 
correcting factor h^, as well as all its building blocks, has the structure 

h[i = i+o{x). 

The reader will find in Ref. [57] the definitions of the quantities entering 
the "Newtonian" waveform h^ e \ as well as the precise definition of the 
effective source factor S^g , which constitutes the first factor in the PN- 
correcting factor h^. Let us only note here that the definition of S^j 
makes use of EOB-dcfincd quantities. For instance, for even-parity waves 
(e = 0) Sgg is defined as the /^-scaled effective energy £ c g/fic 2 . [We use 
the "J- factorization" definition of S^ s when e = 1, i.e. for odd parity 
waves.] 

The second building block in the factorized decomposition is the "tail 
factor" Tg m (introduced in Refs. [M][BB]). As mentioned above, Ti m is a 



resummed version of an infinite number of "leading logarithms" entering 
the transfer function between the near-zone multipolar wave and the far- 
zone one, due to tail effects linked to its propagation in a Schwarzschild 
background of mass Madm = -Heob • ^s explicit expression reads 

_ r(l + 1 - 2lfc) % 2i^lo K (2fcr ) (o 7 \ 

where ro = 2GM/y/e and fc = Gi/lfoBrnfi and fc = mfl. Note that fc 
differs from k by a rescaling involving the rea/ (rather than the effective) 
EOB Hamiltonian, computed at this stage along the sequence of circular 
orbits. 

The tail factor Ti rn is a complex number which already takes into 
account some of the dephasing of the partial waves as they propagate out 
from the near zone to infinity. However, as the tail factor only takes into 
account the leading logarithms, one needs to correct it by a complementary 
dephasing term, e lStm , linked to subleading logarithms and other effects. 
This subleading phase correction can be computed as being the phase Si m 
of the complex ratio between the PN-expanded ftM and the above defined 
source and tail factors. In the comparable-mass case [y ^ 0), the 3PN 822 
phase correction to the leading quadrupolar wave was originally computed 
in Ref. [55J (see also Ref. [53] for the v = limit). Full results for the 
subleading partial waves to the highest possible PN-accuracy by starting 
from the currently known 3PN-accurate ^-dependent waveform [31] have 
been obtained in [57] . For higher-order test- mass (v — > 0) contributions, 
see [53] [52] ■ For extensions of the (non spinning) factorized waveform of 
[57] see [1021 [T0D1 HI5] . 

The last factor in the multiplicative decomposition of the multipolar 
waveform can be computed as being the modulus fg m of the complex 
ratio between the PN-expanded h^l and the above defined source and 
tail factors. In the comparable mass case (y 7^ 0), the /2a modulus 
correction to the leading quadrupolar wave was computed in Ref. |66j 
(see also Ref. [51] for the v = limit). For the subleading partial waves, 
Ref. [57] explicitly computed the other /g m 's to the highest possible PN- 
accuracy by starting from the currently known 3PN-accurate independent 
waveform [31] . In addition, as originally proposed in Ref. [66] . to reach 
greater accuracy the /f m (x;i/)'s extracted from the 3PN-accurate v 7^ 
results are completed by adding higher order contributions coming from 
the v = results |113[ 1114] . In the particular /22 case discussed in [66] , 
this amounted to adding 4PN and 5PN v = terms. This "hybridization" 
procedure was then systematically pursued for all the other multipoles, 



using the 5.5PN accurate calculation of the multipolar decomposition of 
the gravitational wave energy flux of Refs. [1131 1114] . 

The decomposition of the total PN-correction factor into several 
factors is in itself a resummation procedure which already improves the 
convergence of the PN series one has to deal with: indeed, one can see that 
the coefficients entering increasing powers of x in the PN expansion of the 
fern's tend to be systematically smaller than the coefficients appearing in 
the usual PN expansion of h!fl.- The reason for this is essentially twofold: 
(i) the factorization of T(. m has absorbed powers of mir which contributed 
to make large coefficients in h^, and (ii) the factorization of either H c g 
or j has (in the v = case) removed the presence of an inverse square-root 
singularity located at x = 1/3 which caused the coefficient of x n in any 
PN-expanded quantity to grow as 3 n as n — » oo. 

To further improve the convergence of the waveform several 
resummations of the factor ft m {x) = 1 + c[ m x + c| m x 2 + . . . have been 
suggested. First, Refs. [Ml EE] proposed to further resum the fiz{x) 
function via a Pade (3,2) approximant, P| {f22{x\ v)}, so as to improve 
its behavior in the strong-field-fast-motion regime. Such a resummation 
gave an excellent agreement with numerically computed waveforms, near 
the end of the inspiral and during the beginning of the plunge, for different 
mass ratios [M] E2] EI] ■ As we were mentioning above, a new route for 
resumming fg m was explored in Ref. [57 . It is based on replacing fi m by 
its l-ih root, say 

Pl m {x;v) = [f lm {x;v)]^ i . (38) 
The basic motivation for replacing fg m by pi m is the following: the 
leading "Newtonian-level" contribution to the waveform contains a 
factor w^r^ arm w e where /"harm is the harmonic radial coordinate used in 
the MPM formalism [22, 56]. When computing the PN expansion of this 
factor one has to insert the PN expansion of the (dimensionless) harmonic 
radial coordinate rharm, ''harm = a; _1 (l + C\X + 0(x 2 )), as a function of 
the gauge-independent frequency parameter x. The PN re-expansion of 
[niarm(aO] £ then generates terms of the type x~ l (l + lc%x + ....). This is 
one (though not the only one) of the origins of 1PN corrections in hi m and 
fi m whose coefficients grow linearly with I. The study of [57] has pointed 
out that these ^-growing terms are problematic for the accuracy of the 
PN-expansions. The replacement of fg m by pg m is a cure for this problem. 

Several studies, both in the test- mass limit, v — > (see Fig. 1 in [57] ) 
and in the comparable-mass case (see notably Fig. 4 in EZ]), have shown 
that the resummed factorized (inspiral) EOB waveforms defined above 
provided remarkably accurate analytical approximations to the "exact" 
inspiral waveforms computed by numerical simulations. These resummed 



multipolar EOB waveforms are much closer (especially during late inspiral) 
to the exact ones than the standard PN-expanded waveforms given by 
Eq. (|35p with a PN-correction factor of the usual "Taylor-expanded" form 

= 1 + + 4f 2 x^ 2 + 4 m x 2 + ... 

See Fig. 1 in [ST], and slide 29 in my (June 2012) Prague presentation. 

Finally, one uses the newly resummed multipolar waveforms (|34l) to 
define a resummation of the radiation reaction force T v as 

F V = -±F«™ >, (39) 

where the (instantaneous, circular) GW flux F^ max ^ is defined as 

F(£max) = i^EEw 2 w 2 (40) 

1=1 m=l 

Summarizing: Eqs. (fM)) and (f3"9")l . (|40l) define resummed EOB versions 
of the waveform hi m , and of the radiation reaction during inspiral. 
A crucial point is that these resummed expressions are parameter-free. 
Given some current approximation to the conservative EOB dynamics 
(i.e. some expressions for the A,D,Q potentials) they complete the EOB 
formalism by giving explicit predictions for the radiation reaction (thereby 
completing the dynamics, see Eqs. (|2T)j) - (j2"9"j) ). and for the emitted inspiral 
waveform. 

4. EOB description of the merger of binary black holes and of 
the ringdown of the final black hole 

Up to now we have reviewed how the EOB formalism, starting only 
from analytical information obtained from PN theory, and adding extra 
resummation requirements (both for the EOB conservative potentials 
A, Eq. ([21]), an d D, Eq. ([25]). and for the waveform, Eq. ([3"3)l. and 
its associated radiation reaction force, Eqs. (f3"9"]) . (|4"U|) ) makes specific 
predictions, both for the motion and the radiation of binary black holes. 
The analytical calculations underlying such an EOB description arc 
essentially based on skeletonizing the two black holes as two, sufficiently 
separated point masses, and therefore seem unable to describe the merger 
of the two black holes, and the subsequent ringdown of the final, single 
black hole formed during the merger. However, as early as 2000 [H], the 
EOB formalism went one step further and proposed a specific strategy for 
describing the complete waveform emitted during the entire coalescence 



process, covering inspiral, merger and ringdown. This EOB proposal is 
somewhat crude. However, the predictions it has made (years before NR 
simulations could accurately describe the late inspiral and merger of binary 
black holes) have been broadly confirmed by subsequent NR simulations. 
[See the Introduction for a list of EOB predictions.] The original EOB 
proposal (which was motivated partly by the closeness between the 2PN- 
accurate effective metric g^, [30] and the Schwarzschild metric, and by the 
results of Refs. [77] and [109] ) consists of: 

(i) defining, within EOB theory, the instant of (effective) "merger" of 
the two black holes as the (dynamical) EOB time t m where the orbital 
frequency Q(t) reaches its maximum; 

(ii) describing (for t < t m ) the inspiral-plus-plunge (or simply 
insplunge) waveform, /i lnsplunge (i), by using the inspiral EOB dynamics 
and waveform reviewed in the previous Section; and 

(iii) describing (for t > t m ) the merger-plus-ringdown waveform as a 
superposition of several quasi-normal-mode (QNM) complex frequencies of 
a final Kerr black hole (of mass Mf and spin parameter a / , self-consistency 
estimated within the EOB formalism), say 

im) ^ Sd ° Wn ( t ) = E C N e-^-^ , (41) 

with = ajy + ioJN, and where the label N refers to indices (I, l\ m, n), 
with (£, to) being the Schwarzschild-background multipolarity of the 
considered (metric) waveform hg m , with n = 0, 1, 2 . . . being the 'overtone 
number' of the considered Kerr-background Quasi-Normal-Mode, and £' 
the degree of its associated spheroidal harmonics Se' m (aa,9); 

(iv) determining the excitation coefficients of the QNM's in Eq. (l4"Tj) 
by using a simplified representation of the transition between plunge 
and ring-down obtained by smoothly matching (following Ref. |64| ) . on 
a (2p + l)-toothed "comb" (t m ~ p6, . . . ,t m - S, t m , t m + 5, . . . , t m + pS) 
centered around the merger (and matching) time t m , the inspiral-plus- 
plunge waveform to the above ring-down waveform. 

Finally, one defines a complete, quasi-analytical EOB waveform 
(covering the full process from inspiral to ring-down) as: 

hf° B (t) = 9(t m - t) /C PlUnSe W + 9(t - t m ) /C^W . (42) 

where 9(t) denotes Heaviside's step function. The final result is a waveform 
that essentially depends only on the choice of a resummed EOB A(u) 
potential, and, less importantly, on the choice of resummation of the main 
waveform amplitude factor / 2 2 = (P22) 2 - 



We have emphasized here that the EOB formalism is able, in principle, 
starting only from the best currently known analytical information, to 
predict the full waveform emitted by coalescing binary black holes. 
The early comparisons between 3PN-accurate EOB predicted waveform^] 
and NR-computed waveforms showed a satisfactory agreement between 
the two, within the (then relatively large) NR uncertainties [SHJ 1103] . 
Moreover, as we shall discuss below, it has been recently shown that the 
currently known Pade-resummed 3PN-accurate A(u) potential is able, as 
is, to describe with remarkable accuracy several aspects of the dynamics 
of coalescing binary black holes, |96| [73] . 

On the other hand, when NR started delivering high-accuracy 
waveforms, it became clear that the 3PN-level analytical knowledge 
incorporated in EOB theory was not accurate enough for providing 
waveforms agreeing with NR ones within the high-accuracy needed for 
detection, and data analysis of upcoming GW signals. [See, e.g., the 
discussion in Section II of Ref. |100j .] At that point, one made use of the 
natural flexibility of the EOB formalism. Indeed, as already emphasized 
in early EOB work 03 [S3], we know from the analytical point of view 
that there are (yet uncalculated) further terms in the it-expansions of the 
EOB potentials A(u), D(u), . . . (and in the x-expansion of the waveform), 
so that these terms can be introduced either as "free parameter (s) in 
constructing a bank of templates, and [one should] wait until" GW 
observations determine their value(s) [47] , or as " fitting parameters and 
adjusted so as to reproduce other information one has about the exact 
results" (to quote Ref. [S3])- For instance, modulo logarithmic corrections 
that will be further discussed below, the Taylor expansion in powers of u 
of the main EOB potential A(u) reads 

A Taylor (u; v) = 1 - 2u + a 3 (v)u 3 +a 4 (v)u 4 + a 5 (v)u 5 +a 6 (v)u 6 + ... 

where the 2PN and 3PN coefficients a^(y) — 2v and a^v) = a±v are 
known, but where the 4PN, 5PN,. . . coefficients, a^{v) : clq{v) , . . . have not 
yet been calculated (see, however, below). A first attempt was made in 
[54] to use numerical data (on circular orbits of corotating black holes) to 
fit for the value of a (single, effective) 4PN parameter of the simple form 
Q-siy) = a 5 l> entering a Pade-resummed 4PN-level A potential, i.e. 

A\{u- a h ,v) = P\ [A 3PN (u) + isa 5 u 5 ] . (43) 

This strategy was pursued in Ref. |42[ 166] and many subsequent works. 
It was pointed out in Ref. 67 that the introduction of a further 5PN 

8 The new, resummed EOB waveform discussed above was not available at the 
time, so that these comparisons employed the coarser "Newtonian-level" EOB 
waveform h^'^(x). 



coefficient a§{v) — clqv, entering a Pade-resummed 5PN-level A potential, 
i.e. 

A\{u\ a 5 , o 6) v) = F5 1 [A 3PN {u) + ua b u 5 + va 6 u 6 } , (44) 

helped in having a closer agreement with accurate NR waveforms. 

In addition, Refs. [53J introduced another type of flexibility 
parameters of the EOB formalism: the non quasi-circular (NQC) 
parameters accounting for uncalculated modifications of the quasi-circular 
inspiral waveform presented above, linked to deviations from an adiabatic 
quasi-circular motion. These NQC parameters are of various types, and 
subsequent works [72] [7TJ [57] 23] [T3] [HI 1100] have explored several ways of 
introducing them. They enter the EOB waveform in two separate ways. 
First, through an explicit, additional complex factor multiplying he m , e.g. 

/L QC = (1 + 4 m ni + 4 m n 2 ) exp[z(4 m n 3 + 4 m n 4 )} 

where the Uj's are dynamical functions that vanish in the quasi-circular 
limit (with n\, 7J2 being time-even, and 71.3, 724 time-odd). For instance, one 
usually takes n\ = (p r ,/r£l) 2 . Second, through the (discrete) choice of the 
argument used during the plunge to replace the variable x of the quasi- 
circular inspiral argument: e.g. either xq = (GMfi) 2 / 3 , or (following [53] ) 
x v = — (r u Q) 2 where v v = flr u , and r u = r[i)j{r,p ip )] 1 /' i is a modified 
EOB radius, with ip being defined as 



l + 2v \JA(r) l + '-f 




(45) 



For a given value of the symmetric mass ratio, and given values of the 
A- flexibility parameters 05(f), a§{v) one can determine the values of the 
NQC parameters af m 's from accurate NR simulations of binary black hole 
coalescence (with mass ratio v) by imposing, say, that the complex EOB 
waveform hf° B (t EOB ; 05, a^; a\ m ) osculates the corresponding NR one 
/i^(i NR ) at their respective instants of "merger", where t^ B eI = t^ B 
was defined above (maximum of ri EOB (i)), while i„^ rger is defined, say, as 
the (retarded) NR time where the modulus |/i^ 2 R WI °f t ne quadrupolar 
waveform reaches its maximum. The order of osculation that one requires 
between hf° B (t) and h™(t) (or, separately, between their moduli and 
their phases or frequencies) depends on the number of NQC parameters 
of™. For instance, of™ and a^" 1 affect only the modulus of hf° B and 
allow one to match both |/if,^ B | and its first time derivative, at merger, to 
their NR counterparts, while a^^a 1 ™ affect only the phase of the EOB 
waveform, and allow one to match the GW frequency uf® B (t) an d its 
first time derivative, at merger, to their NR counterparts. The above 



EOB/NR matching scheme has been developed and declined in various 
versions in Refs. [73 ED EH S3 QU E2 13 EU US] . One has also extracted 
the needed matching data from accurate NR simulations, and provided 
explicit, analytical ^-dependent fitting formulas for them [671 1100) [70] . 

Having so "calibrated" the values of the NQC parameters by extracting 
non-perturbative information from a sample of NR simulations, one can 
then, for any choice of the A-flexibility parameters, compute a full EOB 
waveform (from early inspiral to late ringdown). The comparison of the 
latter NQC-completed EOB waveform to the results of NR simulations is 
discussed in the next Section. 

5. EOB vs NR 

There have been several different types of comparison between EOB and 
NR. For instance, the early work [39] pioneered the comparison between 
a purely analytical EOB waveform (uncalibrated to any NR information) 
and a NR waveform, while the early work |65j compared the predictions for 
the final spin of a coalescing black hole binary made by EOB, completed by 
the knowledge of the energy and angular momentum lost during ringdown 
by an extreme mass ratio binary (computed by the test-mass NR code of 
[74] ) . to comparable- mass NR simulations [86] . Since then, many other 
EOB/NR comparisons have been performed, both in the comparable- 
mass case |103| S3 EH] E3 EJ EZ] H3] , and in the small-mass-ratio case 
|64[ I121| I122| IT2] [9] . Note in this respect that the numerical simulations 
of the GW emission by extreme mass-ratio binaries have provided (and 
still provide) a very useful "laboratory" for learning about the motion 
and radiation of binary systems, and their description within the EOB 
formalism. 

Here we shall discuss only two recent examples of EOB/NR 
comparisons, which illustrate different facets of this comparison. 

5.1. EOB[NR] waveforms vs NR ones 

We explained above how one could complete the EOB formalism by 
calibrating some of the natural EOB flexibility parameters against NR 
data. First, for any given mass ratio v and any given values of the 
A- flexibility parameters 05(1/), ~a§{i>), one can use NR data to uniquely 
determine the NQC flexibility parameters a^s. In other words, we have 
(for a given v) 

a,; = ai[NRdata; a$, as] , 

where we defined a§ and a§ so that 0,5(1/) = a^v^'a^u) = a§v. [We allow 
for some residual ^-dependence in as and a&.] Inserting these values in the 



(analytical) EOB waveform then defines an NR-complctcd EOB waveform 
which still depends on the two unknown flexibility parameters 05 and olq. 

In Ref. |67] the (05, ag)-dependent predictions made by such a NR- 
completed EOB formalism were compared to the high-accuracy waveform 
from an equal-mass binary black hole (y — 1/4) computed by the Caltcch- 
Cornell-CITA group [112] , (and then made available on the web). It 
was found that there is a strong degeneracy between 05 and a§ in the 
sense that there is an excellent EOB-NR agreement for an extended region 
in the (05, ae)-plane. More precisely, the phase difference between the 
EOB (metric) waveform and the Caltech-Cornell-CITA one, considered 
between GW frequencies Mujl — 0.047 and Mwr = 0.31 (i.e., the last 
16 GW cycles before merger), stays smaller than 0.02 radians within 
a long and thin banana-like region in the (as, ag)-plane. This "good 
region" approximately extends between the points (05, ae) — (0,-20) 
and (05, ae) = (—36, +520). As an example (which actually lies on the 
boundary of the "good region"), we shall consider here (following Ref. [69] ) 
the specific values 0,5 = 0, ag = —20 (to which correspond, when v = 1/4, 
a x = -0.036347, a 2 = 1.2468). [Ref. [57] did not make use of the NQC 
phase flexibility; i.e. it took 03 = 04 — 0. In addition, it used n 2 — r/rQ 2 
and introduced a (real) modulus NQC factor ff^ C only for the dominant 
quadrupolar wave 1 = 2 = m.] We henceforth use M as time unit. 
This result relies on the proper comparison between NR and EOB time 
series, which is a delicate subject. In fact, to compare the NR and EOB 
phase time-series </»^ 2 R (^nr) and 0f 2 OB (tEOB) one needs to shift, by additive 
constants, both one of the time variables, and one of the phases. In 
other words, we need to determine r and a such that the "shifted" EOB 
quantities 

*BOB = *EOB + r , cf>'i OB = 0f 2 OB + a (46) 

"best fit" the NR ones. One convenient way to do so is first to "pinch" (i.e. 
constrain to vanish) the EOB/NR phase difference at two different instants 
(corresponding to two different frequencies uj\ and w 2 ). Having so related 
the EOB time and phase variables to the NR ones we can straigthforwardly 
compare the EOB time series to its NR correspondant. In particular, we 
can compute the (shifted) EOB-NR phase difference 

A— ^ 2 OBNR (*nr) = ^ OB (i' EOB ) - ^ 2 R (i NR ). (47) 

Figure E] compare^ (the real part of) the analytical EOB metric 
quadrupolar waveform ^f^ 5 /^ to the corresponding (Caltech-Cornell- 
CITA) NR metric waveform ^J^P/za [Here, ^ 22 denotes the Zerilli- 

9 The two "pinching" frequencies used for this comparison are Mwi = 0.047 
and Mwj = 0.31. 
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Figure 2. This figure illustrates the comparison (made in Refs. [571 
[69] ) between the (NR-completed) EOB waveform (Zerilli-normalized 
quadrupolar (£ = m = 2) metric waveform (|42l) with parameter-free 
radiation reaction ([55]) and with a 5 = 0, a 6 = —20) and one of the 
most accurate numerical relativity waveform (equal-mass case) nowadays 
available [112] . The phase difference between the two is A</> < ±0.01 
radians during the entire inspiral and plunge, which is at the level of the 
numerical error. 



normalized asymptotic quadrupolar waveform, i.e. ^22 = i?/»22/v / 24 
with R = Rc 2 /GM.] This NR metric waveform has been obtained 
by a double time- integration (following the procedure of Rcf. [71]) from 
the original, publicly available, curvature waveform ipj 2 |112j . Such 
a curvature waveform has been extrapolated both in resolution and in 
extraction radius. The agreement between the analytical prediction and 
the NR result is striking, even around the merger. See Fig. [3] which closes 
up on the merger. The vertical line indicates the location of the EOB- 
merger time, i.e., the location of the maximum of the orbital frequency. 

The phasing agreement between the waveforms is excellent over the 
full time span of the simulation (which covers 32 cycles of inspiral and 
about 6 cycles of ringdown) , while the modulus agreement is excellent over 
the full span, apart from two cycles after merger where one can notice a 
difference. More precisely, the phase difference, Acj) — ^metric — ^metric > 
remains remarkably small (~ ±0.02 radians) during the entire inspiral 
and plunge (o_>2 = 0.31 being quite near the merger). By comparison, 
the root-sum of the various numerical errors on the phase (numerical 
truncation, outer boundary, extrapolation to infinity) is about 0.023 
radians during the inspiral |112j . At the merger, and during the ringdown, 




3800 3820 3840 3860 3880 3900 3920 3940 3960 3980 4000 
t 



Figure 3. Close up around merger of the waveforms of Fig. [21 Note the 
excellent agreement between both modulus and phasing also during the 
ringdown phase. 



Acj) takes somewhat larger values (~ ±0.1 radians), but it oscillates 
around zero, so that, on average, it stays very well in phase with the 
NR waveform whose error rises to ±0.05 radians during ringdown. In 
addition, Ref. |67j compared the EOB waveform to accurate numerical 
relativity data (obtained by the Jena group [71]) on the coalescence of 
unequal mass-ratio black-hole binaries. Again, the agreement was good, 
and within the numerical error bars. 

This type of high-accuracy comparison between NR waveforms and 
EOB[NR] ones (where EOB[NR] denotes a EOB formalism which has been 
completed by fitting some EOB-flexibility parameters to NR data) has 
been pursued and extended in Ref. [100] . The latter reference used the 
"improved" EOB formalism of Ref. [57] with some variations (e.g. a third 
modulus NQC coefficient cij, two phase NQC coefficients, the argument 
x Q = (MSI) 2 / 3 in {pl^ lor {x)) e , eight QNM modes) and calibrated it to NR 
simulations of mass ratios q = m^/mi = 1, 2, 3, 4 and 6 performed by the 
Caltech-Cornell-CITA group [37] HZ]- They considered not only the leading 
(£,m) = (2,2) GW mode, but the subleading ones (2, 1), (3,3), (4,4) and 
(5, 5). They found that, for this large range of mass ratios, EOB[NR] (with 
suitably fitted, ^-dependent values of as and qq) was able to describe 



the NR waveforms essentially within the NR errors. This confirms the 
usefulness of the EOB formalism in helping the detection and analysis of 
upcoming GW signals. 

Here, having in view GW observations from ground-based interfero- 
metric detectors we focussed on comparable-mass systems. The EOB for- 
malism has also been compared to NR results in the extreme mass-ratio 
limit !/<l, In particular, Ref. |12] found an excellent agreement between 
the analytical and numerical results. 

5.2. EOB[3PN] dynamics vs NR one 

Let us also mention other types of EOB/NR comparisons. Recently, 
two examples of EOB/NR comparisons have been performed directly 
at the level of the dynamics of a binary black hole, rather than 
at the level of the waveform. Moreover, contrary to the waveform 
comparisons of the previous subsection which involved an NR-completed 
EOB formalism ("EOB[NR]"), the dynamical comparisons we are going 
to discuss involve the purely analytical 3PN-accurate EOB formalism 
("EOB[3PN]"), without any NR-based improvement. 

First, Le Tiec et al. [HS] have extracted from accurate NR simulations 
of slightly eccentric binary black-hole systems (for several mass ratios 
q = mi/m?, between 1/8 and 1) the function relating the periastron- 
advance parameter 



(where A$ is the periastron advance per radial period) to the 
dimensionless averaged angular frequency M£l v (with M = m\ + rri2 as 
above) . Then they compared the NR-estimate of the mass-ratio dependent 
functional relation 

K = K(MQ, V ; v) , 

where v = q/(l + q) 2 , to the predictions of various analytic approximation 
schemes: PN theory, EOB theory and two different ways of using GSF 
theory. Let us only mention here that the prediction from the purely 
analytical EOB[3PN] formalism for K(MQ, v \v) [35] agreed remarkably 
well (essentially within numerical errors) with its NR estimate for all mass 
ratios, while, by contrast, the PN-expanded prediction for K(M^l lp ; v) [50] 
showed a much poorer agreement, especially as q moved away from 1. 

Second, Damour, Nagar, Pollney and Reisswig [73] have recently 
extracted from accurate NR simulations of black-hole binaries (with mass 
ratios q — mj/mi = 1,2 and 3) the gauge-invariant relation between the 
(reduced) binding energy E = (£ tot — M)/fi and the (reduced) angular 
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Figure 4. Comparison (made in [73]) between various analytical estimates 
of the energy-angular momentum functional relation and its numerical- 
relativity estimate (equal- mass case). The standard "Taylor-expanded" 
3PN E(j) curve shows the largest deviation from NR results, especially 
at low j's, while the two (adiabatic and nonadiabatic) 3PN-accurate, non- 
NR- calibrated EOB E(j) curves agree remarkably well with the NR one. 



momentum j — J j (G/iM) of the system. Then they compared the NR- 
estimate of the mass-ratio dependent functional relation 

E = E{j;v) 

to the predictions of various analytic approximation schemes: PN theory 
and various versions of EOB theory (some of these versions were NR- 
completed) . Let us only mention here that the prediction from the purely 
analytical, 3PN-accurate EOB[3PN] for E(j\v) agreed remarkably well 
with its NR estimate (for all mass ratios) essentially down to the merger. 
This is illustrated in Fig. 4 for the q = 1 case. By contrast, the 3PN 
expansion in (powers of 1 /c 2 ) of the function E(J; v) showed a much poorer 
agreement (for all mass ratios). 




6. Other developments 

6.1. EOB with spinning bodies 

We lack space here for discussing the extension of the EOB formalism 
to binary systems made of spinning bodies. Let us start by mentionning 
that the spin-extension of the EOB formalism was initiated in Ref. [47] , 
that the first EOB-based analytical calculation of a complete waveform 
from a spinning binary was performed in Ref. |38j . and that the first 
attempt at calibrating a spinning EOB model to accurate NR simulations 
of spinning (non precessing) black-hole binaries was presented in 101 . 
In addition, several formal aspects related to the inclusion of spins in 
the EOB formalism have been discussed in Refs. [53] HU 13 EH] IE] (see 
references within these papers for PN works dealing with spin effects) 
and a generalization of the factorized multipolar waveform of Ref. [57] to 
spinning, non-precessing binaries has been constructed in Refs. jl02l 1115] , 

6.2. EOB with tidally deformed bodies 

In binary systems comprising neutron stars, rather than black holes, 
the tidal deformation of the neutron star(s) will significantly modify the 
phasing of the emitted gravitational waveform during the late inspiral. As 
GW's from binary neutron stars are expected sources for upcoming ground- 
based GW detectors, it is important to extend the EOB formalism by 
including tidal effects (see [75] and references therein) . This extension has 
been defined in Refs. [55J [TB]. The comparison between this tidal-extended 
EOB and state-of-the-art NR simulations of neutron-star binaries has been 
discussed in Refs. [2] [3J. It appears from these comparisons that the tidal- 
extended EOB formalism is able to describe the motion and radiation of 
neutron-star binaries within NR errors. More accurate simulations will 
be needed to ascertain whether one needs to calibrate some higher-order 
flexibility parameters of the tidal-EOB formalism, or whether the currently 
known analytic accuracy is sufficient. 

6.3. EOB and GSF 

We mentioned in the Introduction that GSF theory has recently opened 
a new source of information on the general relativistic two-body problem. 
Let us briefly mention here that there has been, recently, a quite useful 
transfer of information from GSF theory to EOB theory. The program of 
using GSF-theory to improve EOB-theory was first highlighted in Ref. [15] ■ 
That work pointed to several concrete gauge-invariant calculations (within 
GSF theory) that would provide accurate information about the 0{y) 
contributions to several EOB potentials. More precisely, let us define the 



functions a(u) and d(u) as the f-linear contributions to the EOB potentials 
A(u; v) and D(u; v) = v): 

) = 1 - 2u + va{u) + 0(v 2 ) , 

D{u] v) = (AB)- 1 = 1 + ud(u) + 0{v 2 ) . 

Ref. [48] has shown that a computation of the GSF-induced correction 
to the periastron advance of slightly eccentric orbits would allow one to 
compute the following combination of EOB functions 

p(u) = a(u) + u a'{u) + — u(l — 2u) a"(u) + (1 — 6m) d(u) . 

The GSF-calculation of the EOB function p{u) was then performed in 
Ref. [B] (in the range Q<u<\). 

More recently, a series of works by Le Tiec and collaborators [5511941 HU] 
have (through an indirect route) shown how GSF calculations could be 
used to compute the EOB v- linear a(u) function separately from the 
d(u) one. Ref. [TU] then gave a fitting formula for a(u) over the interval 
< m < | as well as accurate estimates of the coefficients of the Taylor 
expansion of a(u) around u = (corresponding to the knowledge of the 
PN expansion of a(u) to a very high PN order). Very recently, Ackay et 
al. p] succeeded in accurately computing (through GSF theory) the EOB 
a(u) function over the larger interval < u < 4. It was (surprisingly) 
found that a(u) diverges like a(u) sa 0.25(1 — 3u) -1 ^ 2 at the light-ring 
limit u . The meaning for EOB theory of this singular behavior of 

a(u) at the light- ring is discussed in detail in Ref. [TJ. 

6.4- Toward further improvements to EOB 

Let us finally mention some avenues for further progress in EOB theory. 

Logarithmic contributions to the A(u) and D(u) functions have been 
recently computed at the 4PN level [HI El] and even the 5PN one [49lfT0]. 
They have been incorporated in a recent, improved implementation of the 
EOB formalism [70]. 

Two groups have embarked on a calculation of the (full) conservative 
dynamics at the 4PN level [81] [91] . If they succeed, it will be important to 
translate their gauge- dependent results in the gauge-invariant form used 
in EOB theory. [Remember that EOB theory is essentially based on the 
gauge- invariant Delaunay Hamiltonian H(I a ).] 

More generally, let us emphasize that the EOB formalism provides 
a convenient, gauge- invariant way of packaging both the conservative 



dynamics and the multipolar waveform. This EOB packaging has often 
turned out to be very economical. We recommend that authors computing 
high-order PN corrections to either the dynamics or the waveform 
reexpress their results in terms of the EOB building blocks. 

For instance, Jaranowski and Schacfer [91] have recently given a 
partial result at 4PN, expressed in terms of the (gauge-invariant) function 
S(MO„;i/). In terms of this function, the 4PN contribution is a 
polynomial of the fourth degree in v, namely, with x = (MQ^) 2 ^ 3 and 

E{x;v) = -i/ic 2 x(l + e 1PN (is)x + e 2 pN{v)x 2 + e 3PN (is)x 3 

4 , nf,J> - 



+e4pjv(i/; \nx)x + 0(x lnx)) 



they found 



/ , x 3969 2 301 , 77 4 448 , 

e 4PN ( V ;\nx) = - —+c lV +c 2 u +—v +—v\nx, (48) 

where they could not compute the values of the coefficients c\ and c 2 of 
the terms linear and quadratic in is, but only the contributions cubic and 
quartic in is. We wish to point out that their result is re-expressed in a 
more economical (and more informative) way in terms of the basic EOB 
potential A(u; is). Indeed, in terms of the PN expansion, of A(u; v), 

A Taylor (u; is) = l-2u + a 3 (is)u 3 + a 4 (is)u 4: + a 5 (is;\nu)u 5 + a e (is;\nu)u 6 + . . . 

the information contained in the above result can be entirely re-expressed 
in terms of the 4PN- level coefficient 05 (i^; In it). When doing this re- 
expression, one then finds that the information content of Eq. (|4"5|) is that 
the 4PN-level EOB coefficient a5(^;lnu) is no more than quadratic in v, 
i.e. 

64 

a^iy; In u) — (a 5 + — In u)v + a' b v 2 , 
o 

without contributions of degree v z and v A . We recall that similar 
cancellations of higher v n terms were found at lower PN orders in the 
EOB A(u] v) function. Namely, they were found to contain only terms 
linear in is, while ^(v) could a priori have been quadratic in v, and a,i(is) 
could a priori have been cubic in is. The fact that similar remarkable 
cancellations still hold, according to the result of [91] , at the 4PN level, is 
a clear indication that the EOB packaging of information of the dynamics 
in the A(u; nu) potential is quite compact. Indeed, it says that the two 
complicated terms y^r^ 3 + 31104 ^ 4 m the energy function are already 
encoded in the structure of the EOB formalism. Finally, note that 
the full gauge-invariant content of a 4PN computation of the dynamics, 
when interpreted within the EOB formalism, is described by only three 



EOB terms: the coefficient 0,5(1;; hi u) in A(u;v) 1 an anaiogous coefficient 

G?4(V;lnit) in D(u;v), and an additionai contribution to Q(p). 

Regarding the waveform, iet us mention another recent exampie where 
it would have been useful and clarifying to use the EOB packaging. 
Namely, when re-expressing it in terms of the factorized EOB waveform, 
the new content of the recent 3.5PN level computation by Faye, Marsat, 
Blanchet, and Iyer [79] of the PN-expanded quadrupolar waveform /122, 
is entirely contained in an additional 3.5PN-level contribution to the 
supplementary phase, namely S 22 = (30995/1134 v + 962/135 v 2 ) x 7/2 . 
Indeed, the 3.5PN-level contributions to the modulus computed in[7j5] were 
already included in the factorized EOB waveform of Ref. [57] . 

7. Conclusions 

We hope that this brief review has made it clear that: 

1. There is a complementarity between the various current approaches 
to the general relativistic two-body problem: post-Newtonian, Effective 
One Body, gravitational self-force and numerical relativity simulations (of 
both comparable-mass and extreme-mass-ratio systems). 

2. The effective one body formalism offers a convenient framework 
for combining, in a synergetic manner, information coming from the 
other approaches. This formalism seems to constitute an efficient way to 
analytically describe the motion and radiation of circularized^! binaries, 
and to provide accurate gravitational wave templates for detection and 
data analysis. 

3. The general relativistic two-body problem is more lively than ever. 
It provides an example of Poincare's sentence: "II n'y a pas de problemes 
resolus, il y a seulement des problemes plus ou moins resolus." ["There 
are no (definitely) solved problems, there are only more-or-less solved 
problems."] 
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